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Abstract 

The Low Density Parity Check (LDPC) Code decoding algorithm make use of a scaled 
receive signal derived from maximizing the log-likelihood ratio of the received signal. 
The scaling factor (often called the combining ratio) in an AWGN channel is a ratio 
between signal amplitude and noise variance. Accurately estimating this ratio has shown 
as much as 0.6 dB decoding performance gain. This presentation briefly describes three 
methods for estimating the combining ratio: a ‘Pilot-Guided’ estimation method, a 
‘Blind” estimation method, and a Simulation-Based Look-Up table. The ‘Pilot Guided’ 
Estimation method has shown that the maximum likelihood estimates of signal amplitude 
is the mean inner product of the received sequence and the known sequence, the attached 
synchronization marker (ASM) , and signal variance is the difference of the mean of the 
squared received sequence and the square of the signal amplitude. This method has the 
advantage of simplicity at the expense of latency since several frames worth of ASMs. 
The ‘Blind’ estimation method’s maximum likelihood estimator is the average of the 
product of the received signal with the hyperbolic tangent of the product combining ratio 
and the received signal. The root of this equation can be determined by an iterative 
binary search between 0 and 1 after normalizing the received sequence. This method has 
the benefit of requiring one frame of data to estimate the combining ratio which is good 
for faster changing channels compared to the previous method, however it is 
computationally expensive. The final method uses a look-up table based on prior 
simulated results to determine signal amplitude and noise variance. In this method the 
received mean signal strength is controlled to a constant soft decision value. The 
magnitude of the deviation is averaged over a predetermined number of samples. This 
value is referenced in a look up table to determine the combining ratio that prior 
simulation associated with the average magnitude of the deviation. This method is more 
complicated than the ‘Pilot-Guided’ Method due to the gain control circuitry, but does 
not have the real-time computation complexity of the ‘Blind’ Estimation method. Each 
of these methods can be used to provide an accurate estimation of the combining ratio, 
and the final selection of the estimation method depends on other design constraints. 
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Background & Purpose / 


• The LDPC decoder needs a scaled received signal as its input 

• Scaling Factor used is derived from the log likelihood ratio of the received 

signal 2 A 


> Implies a scaling factor of 2a where a = — 

a 1 

> Though referred to as SNR, it is a combining ratio. In digital communications, symbol SNR 
is: £ A 2 

~N~ 0 ~2^ 

• Without a good estimate there is as much as 0.6dB decoding performance loss 
associated with the LDPC decoding algorithm. 
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Purpose eS&3 


• To compare methods of determining the combining ratio 

> Method 1: The Pilot Guided Estimation Method 

❖ Relies on the known Attached Synchronization Markers (ASMs) at the start of every frame 

❖ Simple implementation 

❖ Necessary to have several frames to make an estimation 

❖ Not ideal in for a fast changing channel 

> Method 2: The Blind Estimation Method 

❖ Makes an estimate with 1 frame of normalized unknown data 

❖ Complicated and iterative process to arrive at the estimation 

❖ Use of an entire frame not ideal for fast changing channels 

> Method 3: Simulation & Look-Up Table Method 

❖ Makes an estimate with 1/2 frame of unknown data 

❖ Requires additional gain control to hold the received signal near a 512 

❖ Can be used for a faster changing channel 
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Simulink Block for Pilot-Guided Estimation Method 
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Variance 


Concatenate 


• Based on the Matlab code used in the analysis of the Pilot- 
Guided Estimation Method (Method 1) 

• The mean and variance of the estimator corresponds to the 
results obtained in analysis. 

• Initialization inputs are set in the Matlab workspace 

• The outputs were verified against analysis 


ESCG/EV7/JSC/NASA 


5 




Avionic Systems Analysis 



The Pilot-Guided Estimation Method 'J 


The Pilot-Guided Method uses a known marker, therefore the Probability 
Density function of a received signal in AWGN (y) is : 


P{y | x) = 




exp 


1 iG 


2(7 


Af 


The derived Maximum Likelihood Estimation’s (MLE’s) from the Log- 
Likelihood Ratio are: 


, ] N 

3 = -Zy,x 
N i = 1 


N i = 1 


y 




The estimator (a) is based on the formula: a 
where (•) denotes mean over N symbols \ / N 

2c is usually the detected signal, however since the ASMs are known 
this is substituted with x 

Analysis will vary the number of ASMs to estimate a at various Es/No 

> 1 ASM to 10 ASM by 1 ASM step 

> Es/No varies from -3 dB to 5 dB in 1 dB increments 
Determine the necessary # of ASM to meet the requirement 

> A Sample Space of 1 0 A 5 ASMs was used 

> Area of Interest occurs when Es/No is around -IdB 

❖ Operating point of the LDPC decoder to output a BER 10 A -8 occurs when the input BER = 10 A -1 

❖ This occurs when SNR = -0.8 dB 
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Mean MLEs of Alpha for different EsNo against true values of Alpha 
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Meas alpha @ Es/No= -2 

■ True alpha @ Es/No= -2 

■ Meas alpha @ Es/No= -1 
• True alpha @ Es/No= -1 

Meas alpha @ Es/No= 0 
True alpha @ Es/No= 0 
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2 3 4 5 6 7 8 

Number of Attached Synchronization Markers (ASMs) 


10 


Es/No (dB) 


-2 

1.217522 

-1 

2.216095 

0 

3.217053 

1 

1.114157 

2.115808 

3.115565 

2 

1.077571 

2.079883 

3.077172 

3 

1.064476 

2.062596 

3.061879 

4 

1.049236 

2.049791 

3.051148 

5 

1.044029 

2.045864 

3.043752 

6 

1.039201 

2.039006 

3.041223 

7 

1.035527 

2.036974 

3.034679 

8 

1.033928 

2.03212 

3.031851 

9 

1.029219 

2.030972 

3.030898 

10 

1.0103 

2.0103 

3.0103 

True 

Values 
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es(?G 


Standard Deviation of MLEs of Alpha for different EsNo 



Es/No (dB) 


-2 

-1 

0 


0.848935 

0.829203 

0.808727 

1 

0.605958 

0.593092 

0.577789 

2 

0.497243 

0.487595 

0.474262 

3 

0.435476 

0.421558 

0.412235 

4 

0.391774 

0.379119 

0.371318 

5 

0.358225 

0.347643 

0.338692 

6 

0.330791 

0.322965 

0.315245 

7 

0.310266 

0.302324 

0.294575 

8 

0.292963 

0.284797 

0.277866 

9 

0.278822 

0.271027 

0.26396 

10 


ESCG/EV7/JSC/NASA 


8 


# of AS Markers 
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Deviation of Combining Ratio Using the Pilot-Guided pr/^p 

Estimation Method 


Deviation of Alpha from True Value for EsNo = -1 dB 


I Within 3 sigma 
Within 2 sigma 
| Within 1 sigma 



23456789 
Number of Attached Synchronization Markers (ASMs) 
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Summary of the Pilot-Guided Estimation Method (2 


• It is guaranteed that all estimates within 1 standard deviation of 
the mean estimate (68%) will meet the requirement of being 
within 0.5dB of the true value 

• Some estimates within 2 standard deviations (95%) may NOT 
meet this requirement 
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The Blind Estimation Method 
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Assumes that there is equal probability (in a BPSK scheme), therefore 
the Probability Density function of a received signal is: 


N 1 

P(y) = rij 

i = 1 z 


-V2 


_ e -{yi- A f/ 2 ( r 2 


-{y i+ Af/2c 


n 




;r 


The Blind Estimation Method uses the Log-Likelihood Ratio to find an 
MLE for the signal amplitude, which is: 

~ 1 ^ jAyA 

^ = ^2>. tanh -T 

A i = 1 V & ) 

Using a single frame, step through an algorithm to estimate 

> The iterative algorithm is used to find a numeric approximation 

> Noise power is calculated once a signal amplitude is determined 

Determine the necessary iterations of the algorithm to meet 
requirements 

> A sample space of 10 A 5 frames was used 

> Area of Interest occurs when Es/No is around -IdB 

❖ Operating point of the LDPC decoder to output a BER 10 A -8 occurs when the input 
BER = 10 A -1 

❖ This occurs when SNR = -0.8 dB 
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The Blind Estimation Method Algorithm (/ 



Where: 

❖ A is the estimate of the input signal amplitude 

❖ Y| is the amplitude of the received bit 


❖ a 2 is the noise power 

❖ N is the number of bits in the frame 


• 1) Normalize the received frame of bits such that the average power of a 
received bit is 1 

> This allows noise power to be defined: o- 2 =1 -A 2 mid 

> Limits the range of estimated signal amplitude to be between 0 and 1 

• 2 ) Set A min near 0 and A max =1 -A min . 


• 3) Calculate the midpoint of the expected range and set this as the estimated 
signal amplitude. 

• 4) Calculate noise power from signal amplitude. 

• 5) Evaluate the equation above with these calculated values. 

• 6) If step 5 is greater than 0 set A max to the midpoint, otherwise set A min to the 
midpoint. 

• 7) Return to step 3 for the proscribed number of iterations. 

• 8) Take the midpoint one final time and set that as the expected signal 


Note that A min cannot be 0 as this generates the trivial solution A = 0. 


amplitude (A mjd ). Use this value to calculate: a 
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Method 2: Why Step 6? O 


F(x) vs. x for 20 Trials w/ N=1024, A = 0.783, & EbNo = -3 



F(x) vs. x for 20 Trials w/ N=1 024 , A= 12, & EbNo = -2 
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Simulink Blocks for Blind Estimation Method 
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Constantl 


• Based on the Matlab code used in the analysis of the Blind 
Method (Method 2) 

• The mean and variance of the estimator corresponds to the 
results obtained in analysis. 

• Initialization inputs are set in a Matlab workspace. 

• The outputs were verified against analysis 
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Mean MLEs of Alpha for different EsNo against true values of Alpha 



4 5 6 7 

# of Iterations in Algorithm 
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Standard Deviation of a with Floating Point Output 
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Standard Deviation of Alpha for various EsNo 



Es/No (dB) 


-2 

-1 

0 

0.055095 

0.052422 

0.049486 

2.075299 

0.194406 

0.049486 

1.147862 

0.128529 

0.049486 

0.602455 

0.114091 

0.645361 

0.34819 

0.337665 

0.345355 

0.303567 

0.26505 

0.238655 

0.293566 

0.249758 

0.218892 

0.291016 

0.24579 

0.21299 

0.290417 

0.244789 

0.211533 

0.290258 

0.244584 

0.211125 
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# of Iterations through the algorithm 
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Deviation of Combining Ratio Using the Blind 

Estimation Method 
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Deviation of Alpha from True Value for EsNo = -1 dB 


Within 3 sigma 
Within 2 sigma 
Within 1 sigma 


-0.5 


4 5 6 7 

Number of Iterations 
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Summary of the Blind Estimation Method / 


• It is guaranteed that all estimates within 1 standard deviation of 
the mean estimate (68%) will meet the requirement of being 
within 0.5 dB of the true value 

• At least 9 iterations of the binary search are needed before 
additional iterations are negligible as compared to the variance 
of the data 
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Prior Simulation - Look-up Table Method '2 


Uses prior simulation to determine: 


k* E< 


1-M 


A, 


• For each SNR the expected value is 
determined by simulation and placed in 
the look-up table 


• The deviation around the specified mean, S i 2 
k, is averaged over the same number of 
samples 


• The average (SNR ADD in the diagram) is 
sent to the look-up table to determine the 
combining ratio 


Scaled samples 



20 Oil Accumulate and Dump 


• It is important to control the strength of 
the received signal in order to accurately 
estimate the signal deviation. 


March 30, 2012 
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MATLAB Simulation: 
Combining Ratio Estimate 



• In the AWGN channel block Es/No was applied from -1 .5 dB - 8.5 dB in 
0.5 dB steps 

• Statistics were calculated over 100000 estimates of the Combining 
Ratio, estimated from 1024 soft decision bits. 


March 30, 2012 
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Automatic Gain Control 
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Absolute Value of 
Scaled samples 



• Provides feedback to keep the received signal within the range 
of the estimation 

• This will adjust the gain of the input samples based on the 
deviation from the baseline 512 point 

• The deviation is accumulated and the 11 most significant bits is 
the ‘Soft AGC’ factor 


March 30, 2012 
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Results 


Combining Ratio Bias (dB) 



Standard Deviation of Means from Estimated Mean Combining Ratio (dB) 



• Bias is within 0.15 dB 

• Standard Deviation (STD) is within 0.5 dB 
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Effects of the Number of Symbols used to pc^p 
Calculate the Combining Ratio 


• Bias is close to unchanged 

• The STD decreases with an 
increase in the number of 
symbols used for the 
Combining Ratio estimation 

• From 128 symbols to 256 
symbols = 0.25dB 
improvement in the STD 

• From 1024 symbols to 2048 
symbols there is a O.ldB 
improvement in the STD 




Combining Ratio Bias (dB) 


128 symbols 
256 symbols 
512 symbols 
1024 symbols 
2048 symbols 


3 4 

Es/No (dB) 

Standard Deviation of Combining Ratios (dB) 


3 4 

Es/No (dB) 
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Summary of Prior Simulation - Look-up 

Table Method 


es(?G 


• The Combining Ratio Estimate is within 0.5 dB of its true value for 71 .2% 
of the data at Es/No = -1 dB (Operating Point) 

> The Combining Ratio Estimate is within 0.5 dB of its true value more than for 98% 
of the data at Es/No = 8.5 


March 30, 2012 
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Summary 8S(^G 


• The three methods for estimating the Combining Ratio - the Pilot- 
Guided Estimation Method, the Blind Estimation Method , and a Prior 
Simulation - Look-up Table Method - are within 0.5 dB for one STD of 
the CR at Es/No = -1 dB. 

• The Pilot-Guided Estimation Method was easy to implement; however, 
6 frames of data were required to estimate one STD of the CR within 
0.5 dB of the true CR. 

• The Blind Estimation Method was complicated to implement and one 
frame of data estimated one STD of the CR within 0.5 dB of the true 
CR. 

• The Prior Simulation - Look-up Table Method using 1024 sample 
estimated one STD of the CR within 0.5 dB of the true CR. 



ESCG/EV7/JSC/NASA 




Avionic Systems Analysis 


References 6S(zG 


• Andrews, Kenneth, et al. “Symbol Scaling for LDPC Decoders.” February 
6,2008. 

• Li, Bin, et al. “A Low Bias Algorithm to Estimate Negative SNRs on an AWGN 

Channel.” IEEE Communications Letters, vol. 6, no.11, November 

2002. 

• Gagliardi, Robert M., Thomas, C. Melvil. “PCM Data Reliability Monitoring 
Through Estimation of Signal-to-Noise Ratio.” IEEE Transactions on 
Communication Technology, vol. com-16, no. 3, June 1968. 

• Performance Evaluation and Modeling of ‘SNR’ Estimators (EV7-09-4473) 

• Modeling and Performance Evaluation of the Combining Ratio Estimator in the 
Orion S-Band Transponder (EV6-1 2-4732) 

• Combining Ratio Table Lookup for LDPC Decoder (Tech Memo 12192-TM026) 

• Soft Decision Processing (Tech Memo 12192-TM022) 

• Symbol Scaling for LDPC Decoders: Kenneth Andrews, Sam Dolinar, Dariush Divsalar, 
Jon Hamkins, Vic Vilnrotter 



ESCG/EV7/JSC/NASA 




Avionic Systems Analysis 


es<?e 



Extra Slides 


ESCG/EV7/JSC/NASA 


27 




Avionic Systems Analysis 



Maximum-Likelihood Estimates (MLEs) ^ 


• The MLE maximizes either Likelihood Ratio or the Log-Likelihood Ratio 


• MLE of a ratio is equal to the ratio of MLEs 


a = 


A 


<j 


Method 1 uses a known marker, therefore the Probability Density function of a 
received signal (y) is : 


P(y I x) 


1 




exp 


1 i(y,x-Af 


2 C7~ i = 1 


Method 2 assumes that there is equal probability (in a BPSK coding scheme), 
therefore the Probability Density function of a received signal is: 


N 1 

ny)=n - 2 

i = 1 z 


,2 


1 c -[ yi - A ) 2 1 2 ° 2 | 1 g ~(> ’ i + A ?/ 2 

cr4l 71 aV 2 71 


We find the maximum by setting the derivative to 0 
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Determining Log-Likelihood Ratio for the Pilot-Guide^ 

Method 
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Determining Log Likelihood Eguetion for the Blind a 

Method pt. 1 k ~L m 
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Determining Log Likelihood Eguetion for the Blind a r /17 a 

Method pt. 2 
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